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ABSTRACT 


Optimal maintenance policies are designed to reduce the number of system failures 
and minimize the cost of repair by scheduling planned replacements. In this area the 
problem of updating the maintenance policy using the past maintenance history has not 
been adequately solved. In this thesis we study a sequential estimation procedure in a 
nonparametric setting to estimate the age replacement policy that minimizes long run 
expected maintenance costs. 

This thesis begins with the discussion of the concepts of preventive maintenance, age 
replacement policies, the settings of our simulation model, and a detailed description of 
the sequential estimation procedure. We include examples using actual replacement data 
which demonstrate the usefulness of the sequential estimation procedure. Monte-Carlo 
methods are used to study the behavior of estimated optimal age replacement policy for 
different sample sizes, costs and underlying system life distributions. We also make 
comparison with Frees and Ruppert’s (1985) sequential procedure for estimating optimal 
age replacement policies. These comparisons show that our sequential estimation pro- 
cedure is competitive and for large sample sizes performs better than the Frees and 
Ruppert’s procedure. Finally, we will introduce a graphical method to estimate the op- 


timal age replacement policy. 
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I. INTRODUCTION 


A. BACKGROUND 

Optimal maintenance policies are designed to reduce the number of system failures 
and minimize the cost of repair by scheduling times for replacements. By far, most of 
the research in this area has been from the modeling standpoint. Even in the most basic 
scenario, where the underlying system lifetimes are assumed to be independent and 
identically distributed, the problem of updating the maintenance policy using the past 
maintenance history has not been adequately solved. In this thesis Monte-Carlo meth- 
ods are used to study a particular nonparametric procedure which updates estimates of 
an optimal age replacement policy after each replacement. In particular, we will focus 


on finding conditions under which this sequential procedure does well. 


B. MAINTENANCE POLICIES 

The effectiveness of a working system depends not only on the innate properties 
built into it in the design and production stages, but also on the quality of its operation, 
maintenance and repair. Maintenance policies are designed to reduce the incidence of 
system failure by scheduling maintenance actions. For example, a maintenance policy 
is warranted when the failure of a system during actual operation is costly or dangerous 
(i.e., When the cost C, of unscheduled maintenance due to system failure is more than 
the cost C, of scheduled maintenance before svstem failure). If the system is character- 
ized by a failure rate that increases with age, it may be wise to replace it before it has 
aged too greatly. This is the concept of preventive imaintenance. Broadly speaking, 
preventive maintenance is the total of all service functions aimed at maintaining and 


improving reliability performance characteristics and concerns itself with such activities 




















as the replacement and repair of systems, inspections, testing and checking of working 
parts during their operation. 

In the preventive maintenance model developed here, it will be assumed that the 
maintenance action returns the equipment to the “as good as new” condition, thus 
providing the same services as equipment that has been “replaced”. This assumption 
implies that the times between failure are independent and have the same distributions. 
If this is not the case, the model needs to be modiuied. In addition, the model assumes 
that the scheduled and unscheduled replacement costs C,, C, remain constant over time. 

The two most common replacement policies are the policy based on age (age re- 
placement) and the policy based on time (block replacement). Under age replacement 
a system is always replaced at the time of failure or at age T, whichever occurs first. 
Under a block replacement policy the system is replaced upon failure and at times T , 
2T , 3T,..., etc. By its nature, age replacement is administratively more difficult to 
implement, since the age of the unit must be recorded. On the other hand, block re- 
placement, although simpler to administer since the age of the unit need not be recorded, 
leads to more frequent replacement of relatively new items [Ref. 1: pp. 178-182]. Thus, 
in this thesis we only consider age replacement policies. 

We define an optimal age replacement policy as the age replacement policy which 
yields the smallest long run expected replacement costs per unit time. To determine an 
optimal age replacement policy we require explicit knowledge of the system’s life dis- 
tribution. But in a real world scenario, we may not know the life distribution explicitly, 
thus the optimal age replacement policy needs to be estimated. Estimation when a fixed 
sample of independent and identically distributed (/.i.d.) system life times were availuble 
has been examined in detail [Refs. 2,3,4,5]. To gather such iid. data, experimental sys- 


tems must be left in service until failure. Therefore, the experimenter can not implement 











an age replacement policy to achieve cost savings while collecting data for estimation. 
When the luxury of observing system lifetumes until failure is not available, a more cost 
effective approach is to estimate sequentially. By this we mean that the estimator of the 
optimal age of replacement is updated after each system replacement. In addition, a 
system under observation is subject to a replacement policy that is close to the best es- 
timated policy so far. The nonparametric procedure of this type which we will study is 
described in detail in Chapter 2. The difficulty with such a sequential procedure is that 
implementing an age replacement policy while collecting data results in censored obser- 
vations. Thus, there is a trade-off between the goals of controlling cost and gathering 
data for estimation. This conflict is particularly acute in the nonparametric setting 
Where information about the right tail of the distribution needed to estimate the optimal 
replacement age can not be obtained under heavy censoring. 

This thesis will include small, moderate and large sample studies by Monte-Carlo 
simulation to determine the behavior of the sequential procedure under various condi- 
tions. In Chapter 2, description of the age replacement policy and the sequential esti- 
mation procedure are given. Included in tlis chapier is an example of the sequential 
estimation procedure applied to actual data. Chapter 3 establishes the setting for the 
simulation study, with a quick overview of the results of that study. A detailed analysis 
of the simulation results is given in Chapter 4, In that chapter, we also present results 
comparing our sequential estimation procedure with another sequential estimation pro- 
cedure. In Chapter 5, we describe a graphical method for estimating the optimal age 


replacement policy. Conclusions and recommendations are given in Chapter 6. 


be 











Il. THE SEQUENTIAL ESTIMATION PROCEDURE 


A. THE OPTIMAL REPLACEMENT AGE 

Let (X,} be an independent and identically distributed (.i.d.) sequence of positive 
random variables with distribution function F. The sequence .¥), ¥,.... represents the 
sequence of system lifetimes that would be observable if the systems were replaced at 
failure. [t is intuitively clear that it is not advantageous to use an age replacement policy 
for a system whose lifetime has a decreasing failure rate or constant failure rate since this 
would result in a replacement by a new system whose properties are worse than or just 
the same as the original one. Let C, be the cost of an unscheduled replacement (at sys- 
tem failure) and C, be the cost of a scheduled replacement (before system failure), where 
C,>C, Under an age replacement policy with scheduled replacement at age , standard 
results from renewal theory [Ref. 6: p. 87] can be used to show that the long run ex- 


pected cost per unit time is given by 


C, x F(t) + C, x S(a) 


R(t) = : 
| S(x) dx 
0 


where S(x) = 1 — F(x) is the survival function. The numerator is the expected cost of 
one replacement under the age replacement policy, and the denominator is the expected 
time between replacement. Under some fairly general conditions, there exists unique and 
finite time @* where R(s) attains a global minimum [Ref. 7: pp. 161-168]. For example, 
a sufficient condition for the existence of @* is that F have failure rate A(x) that strictly 


increases to infinity. We will assume through out the sequel that ™ exists. 














When such a $* exists, it is the optimal age of replacement, in the sense that a 
policy with age replacement at @* will have minimum long run expected costs per unit 


time. When $* is unknown then it needs to be estimated from system lifetime data. 


B. THE SEQUENTIAL ESTIMATION PROCEDURE 

Let {@*} be the sequence of estimators of @* where the estimator }* is based on 
data from the first # replacements. If our goal was just to estimate @*, then we would 
construct @* based on observing ” complete system lifetimes. In this case all replace- 
ments Would be unscheduled replacements. at system failure, with cost C,. We wish, 
however, to control replacement costs by implementing age replacement policies while 
collecting data to estimate @*. There is clearly a trade-off between gathering data for 
estimation (i.e, observing systera. as long as possible) and controlling costs by imple- 
menting the age replacement policy with the best estimated scheduled replacement age 
so far. Therefore, we compromise and chceose the scheduled replacement age for the 
ath system to be €, where 6%,<&,<oo. Intuitively, it scems reasonable to take 
é,= o*, +a, where {a,} is a sequence of. stants that decrease to zero as 00. How 
to choose the sequence {a,} is an open question that we study in detail in the Monte- 
Carlo simulations that follow. 

Assuming that the {a} sequence and ¢,, the replacement age for the first system 
under observation, are given, the procedure to compute the estimators { $* } developed 


in [Ref. 8] follows: 


1, Determine ith system’s scheduled repiacement time by 


= diy + y: 














2. At the ith replacement observe the system lifetime X, or €, whichever comes first. 
Let Z, = min(X,, ¢,) and 6, = I(X, < €,) where J(A) is indicator function of the set A. 
In other words, if system is replaced before failure, then the replacement time 
Z,= é, and 6, = 1, otherwise Z, = X, and 6, = 0. 


3. Sort Z,, Z,,..,Z, into an ascending sequence Zi, Zea) very Ziq With Say) Seyy voy Oey 
ordered according to the ordering Of Z.), Ze), «5 Zea 


4, Use the data (Z)), dq), (Zayy 5a), «+1 (Ziy 5g) to estimate S=1—F, The estimator 
S of S is defined as 


Note that although S is formally identical to the Kaplan-Meier estimator 
[Ref. 9: pp. 34-35] based on randomly right censored data, it’s properties are sub- 
stantially different because the underlying data pairs (Z,, 6,) , (Z,,6;), «.. , (Z,, 5) 
are not independent. 


5. Use the estimator in step 4 to estimate R(1) as 


A A 
A C,x FO} + Cx S(t 
R(t) = 1X O 2X (t) ; 
| S(x) dx 
0 


6. The minimum of the estimator R(s) does not exist. In addition if by = 0, ie, if the 
largest replacement age observed so far occurred before system failure, ther 
R(t) +0 asn->co. To take care of these two problems, after i replacements, We 
estimate @* by o*, where $* attains the infimum of R(t) overo<t< Zp, In fact 
@* is easy to calculate. It can be shown [Ref. 8], that 6 is equal to the Z, among 


1<j<isuch that 


R(Z) = min { R(Zq), R(Zq) 5» R(ZG) }- 


The procedure is then repeated. 











This sequential estimation scheme means that after (i— 1) replacements the exper- 
imenter has available a scheduled replacement age €, The replacement cost for the ith 
system is C, if X¥,<€, otherwise the replacement cost is C, Thus with this sequential 
estimation procedure the actual total replacement cost for the first 7 systems that are 


observed is 


= 


i 


Ch = 


{Cy x 6, + C, x (1-4;)}, 
F 


1 


and the total operating time for the ” systems is 


n 


n 
m= dD, min(X, &) =) 

i=1 = 
Let e>0, and a,>e fori = 1, 2,.... Then in (Ref. 8] it is shown, that as 1 00, 6* 
converges to #* with probability 1.0, and the actual cost per unit time < converges 
with probability 1.0 to R(f*). Since e>0, can be arbitrarily small, for practical pur- 
poses we can choose a sequence {a,} such that 4,90. It is clear that choice of {a,} can 
have considerable impact on how well ¢* estimates ¢* and on how large the actual 


costs oi are. An important feature of this thesis is to use simulation to provide practical 
n 


guidelines for choice of {a,}. 


C. EXAMPLE: REPLACEMENT COST ANALYSIS OF THE TRACTOR AT 
ENGINE FAILURE 

To show how the sequential procedure works, we consider the data 
([Ref. 10: p. 560], summarized in Table 1, which includes the age at engine failure of 22 


brand new D9G-66A Caterpillar tractors and the calendar date of the failure. 

















Assume the unscheduled replacement cost C, = $200.00 and scheduled replacement 
cost C, = $100.00. If engines are renlaced only at failure, with no age replacement 
policy (ie., €,= 00 for i= 1, 2,...., 22), then at each replace:. Z,= X, and F becomes 
the empirical distribution based on the system lifetimes observed so far. Table 2 on 


page 9 gives the estimates of @* calculated after each replacement. 


Table 1. THE AGE OF TRACTOR FAILURE 


‘ : Age in hours when 
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From Table 2, we have 





total replacement cost = 22 x 200.00 = 4400.00, 


22 


total operating hours = > Z, = 135299, and 
f=} 


4400.00 = 0.03252 dollars 


cost per unit time 


135299 hour 

















We now use the sequential estimation procedure described in the last section on the 


same set of data, with the same costs C, = 200.00, C, = 100.00. For this example, let 
500 
i- 





€, =o, so that first replacement will occur at system failure and let a, = hours 


- | 
fori = 2,3,...., 22. The sequence {a,} used in this example was chosen for convenience. 


No special attempt was made to find the {a,} sequence that gave the best results. Ap- 
plying the sequential estimation procedure introduced in last section, the sequential es- 


timates of @* and the data are summarized in Table 3. 


Table 3. ESTIMATION WITH AGE REPLACEMENT POLICY 
a i ee ee 
OR A ae A IT | ACT TC 
2 [300.00 | se61.00 | sos | 1 | s78.00 | sor 
$411.09 
332.67 | 67s | 0 | sx767 | sto 
7 [#333 | seaas | 3556 | 0 | sagas | stor | 
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[ao [362s [2 T2008 


25. 0 5110.00 ee 1 4815.00 5085 


ae [aa s1061 [oss [0 [sos] 085 
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This procedure yields 





total replacement cost = », (200.00 x 6, + 109.00 x (1-6,)) = 3000.00 , 
i=] 
total operating hours = » Z, = 107395.48, and 
cost per unit time = STOEET = 0.02793 collars. : 
At each updating, both procedures give identical estimates ¢*, / = 1, 2, ...., 22. But 


the final cost per unit time when no age replacement policy was implemented is larger. 

By varying unscheduled maintenance cost, it can be seen in Table 4 that as the ratio of 

om 

G_ increases, the benefit of applying an age replacement policy while gathering data 
2 


from estimation becomes greater. 


Table 4. COMPARISON OF MAINTENANCE COSTS 































Average 


























Average 



















Sched- ; Estimated | | . cost with 

uled re- se babi ae optimal eae es Failure 
replace- | placement] ,. ). : replace- a Replace- 
ment cost DOUSUIIE: (HANCE-COSE ment time placement ment Pol- 


cost Policy icy 


| 300 | 100 | 106358.44 3800 | 4813 03573 | 0.04878 
[700 | 100 oes 
9 | Ypres | 260 [oes [oso 
wo frases || 50 [oor 














11 








fil. SIMULATION SETTING 


A. UNDERLYING LIFE DISTRIBUTION 

Guidelines for choosing a sequence {a,} for estimation will depend on three factors, 
the underlying distribution F, the costs C,, C, and the sample size V where we intend to 
stop sampling. The two factors that we assume are known, are the sample size .V and 
the costs C, and C,. In much of the simulation, we concentrate on fixed costs C, = 5.0 
and C, = 1.0. Other costs are also considered, but in much less detail. We choose 
sample sizes V = 10, 50 and 250 to reflect small, moderate and large sample sizes. The 
third factor is the underlying system life distribution which in general will be unknown. 
For simulation, we choose this distribution to be Weibull with shape parameter « and 


scale parameter /, where the density is given by 
fe) = od (A027 eo fort>0. 
The Weibull distribution has failure rate 
A(t) = aA (An. 


When « > 1.0, the failure rate is strictly increasing to infinity. Thus, for Weibull dis- 
tributions, with a > 1.0, a unique and finite optimal replacement age $* exists. The ten 
different Weibull distributions used in the simulation have « values 1.1, 1.2, ...., 1.9, 2.0. 
This selection of « values gives us a range of distributions which become more like the 
exponential distribution as a decreases from 2.0 to 1.1. To make fair comparisons be- 
tween Weibull distributions, the scale parameter 4 was chosen so that the expected sys- 


tem lifetime E(X,) = 2.0. See Figures i and 2, for plots of the Weibull densities and 


corresponding faiiure rates. 
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Figure 1. The Weibull density function f(t) with E(X,) = 2.0 
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, Figure 2. The failure rate of Weibull distribution with E(X,) = 2.0 











B. OPTIMAL REPLACEMENT TIME 
When the distribution of system lifetime is Weibull, the expected long run average 
cost per unit time under a simple age replacement policy with scheduled replacement at 


age t, is given by 


C,(1 - gr ac, ga (ary 


f 2 
| eT AY dy 
0 


Ri) = 


See Figure 3 for a plot of R(t) when the underlying life distribution is Weibull wicn shape 
parameter ¢ varying from 1.1 to 2.0. For each curve on Figure 3 the optimal replace- 
ment time @* can be located on the x-axis at the minimum point. 

Table 5 on page 16 gives the optimal replacement times $* for different values of 
the shape parameter « when the unscheduled replacement cost C, = 5.0, and the 
scheduled replacement cost C, = 1.0, The optimal replacement times #* in Table 5 are 
computed from simulation results.1 Each @* was generated using 60,000 pseudo ran- 
dom Weibull lifetimes. Included in Table 5, is the probability that a system will be re- 


placed before failure under the optimal age replacement policy, 


* vs 
P(X,< 6%) = 1-e7 7%), 


From Table 5 on page 16, we observe that the optimal replacement time * in- 
creases as the shape parameter a decreases from 2.0. For values of « close to 1.0, the 
performance of a new system will not differ greatly from the old one which is still in use. 
In this case, very little is gained by replacing the system before failure at the higher cost 
C, The larger values of }*, insure that a small percentage of replacements will be made 


1 Since the optimal replacement time ¢* comes from a simulation result, it varies slightly 
with the number of pseudo random variables used and the seed numbers used to generate them. 
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Figure 3. Long run expected average cost curves with E(X,) = 2.0 


before failure, which is what we desire if the system’s life distribution is close to expo- 


nential. 


C. OVERVIEW OF SIMULATION RESULTS 
A wide range of {a,} sequences have be on selected for study. Their functional forms 
are given in Table 6 on page 18. Recall that &, (i.e., the scheduled replacement time or 


scheduled censoring time for the ith system lifetime) is 


f= oh, + 4%, b= 2; 3h 3 
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Table 5. ESTIMATED OPTIMAL REPLACEMENT TIMES OF WEIBULL 
MODEL WITH E(X;) = 2.0 


Long run expected 
a a a a 
[19 eamiess [rises | veroons | oasis 
ed 
Tn [eae | aoe [20a | oon 


We have arbitrarily selected €, = 1.0 and let a, = 0.0 for all {a,} sequences. For each 














a , sample size V and sequence {a,} , we have simulated the results of sequentially esti- 
mating @*. Each simulation is based on generating 1000 sequences of system lifetimes. 

Let Ry, j = 1, 2, ..., 1000 be the actual replacement cost per unit time for the first 
N replacements of the jth repetition of a simulation. For each simulation, the perform- 


ance of the sequential estimation procedure is first evaluated by computing 


—  -Y (Ry— RO")? 
MSE = 2, Sea 117 














where MSE is the average squared difference of the actual replacement cost per unit time 
from the minimum long run expected replacement cost per unit time. Because the ulti- 
mate goal of estimating @* is to reduce costs, we have chosen to evaluate performance 
by comparing at actual costs per unit time to R(d*), rather than looking at the mean 
square error of the estimators of @*. Tables 7, 8 and 9 summanie the simulation results 
of the {a,} sequences that performed best for different values of « and sample sizes .\, 
with fixed costs C, = 5.0, C,=1.0. The criteria for selecting the best performing {a,} 
sequence is to choose the sequence which gives the smallest MSE. In cases where {a,} 
sequences had nearly equal MSE, we select the sequence with the lowest actual cost per 


unit time averaged over the 1000 repetitions, 


1000 
R. a Ry 
y= Thon, . 
V™ £, 1000 
Jal 
We have also plotted the best performing {a.} sequences using the three sample sizes 
of 10, 50 and 250 with various values of the shape parameter « in Figure 4 on page 22. 
From these plots and the results frora Tables 7, 8 and 9, we conclude that in general 
when « decreases from 2.0 (ie., the underlying life distribution is becomes more expo- 
nential), the best performing {a,} sequence tends to have larger values. In the next 
chapter we provide a detailed analysis of these simulation results, that reveals how the 
best performing {a,} sequences are effected by sample size, underlying distribution and 


costs. 
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Table 6. FUNCTIONAL FORM OF {a,} SEQUENCES 


Functional form 


Ronee! ieee 
AA cod 


ee, ee 
(A, +4, x (1 —1)) 


] 
A, + (A, x (I — 1) )% 
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Teble 7, BEST SIMULATION RESULTS FOR SMALL SAMPLE SIZE (N= 10) 


Esti- 
MSE of 
o* 

time or 


2.0 pee tee ele ee 2.01557 | 0.36262 | 1.26406 0.27861 
1.725+0.06 x (i— 1) 

: Seen | eae 2.06978 | 0.38469 0.28013 
1.725 +0.06 x (i— 1) 

Seat: Beer eee 2.13750 | 0.41103 | 1.27617 0.32854 
1.725+0.01 x (i— 1) 

ee ets ee =. 2.21857 | 0.44410 | 1.33542 0.38985 
1.400 + 0.03 x (f= 1) 

2.29646 | 0.49618 | 1.36779 0.41863 

0.54815 | 1.39615 0.45544 
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Table 8. BEST SIMULATION RESULTS FOR MODERATE SAMPLE SIZE 


(N= 50) 
Esti- 
mated 
optimal MSE of 
replace- o* 
ment 
time o* 
1.93167 | 0.08290 | 1.26858 0.17087 
1.99418 | 0.08626 | 1.30783 0.19786 
——___! | 2.05713 | 0.08936 | 1.33833 0.22572 
(9.00+ 0.03 x (i= 1))?5 
17 eA tos = So NS ee 2.13354 | 0.09605 | 1.38420 0.25644 
2.225+0.03 x (i— 1) 
i Ns ce Il SP IOTSS ‘ONO 206 0.34657 
(8.00 +0.03 x (i — 1) 







































Actual MSE of 
aie average | actual 
{a,} functional form cost average 
Ry cost 







! 
2.225+0.06 x (i— 1) 














1 
2.525 +0.03 x (i— 1) 





1s | 2.30572 | 0.10254 | 1.59010 | 0.41487 
1.400 +0.03 x (i— 1) 
——___!___ | 2.38234 | 0.11150 | 1.74888 
1.400 +0.03 x (7— 1) 


















0.55554 
1 


1 
1.400 +0.03 x G— 1) 


2.47548 | 0.12823 | 1.91380 0.8325 
2.56903 | 0.15299 | 2.16169 3.34231 
2.68135 | 0.20829 | 2.35770 | 113.78398 
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I 
1.400+ 0.03 x (/— 1) 


1 
1.400 + 0.03 x (i— 1) 






























Table 9. BEST SIMULATION RESULTS FOR LARGE SAMPLE SIZE 
(N= 250) 












































Actual MSE of 
Shape eu 
param- | {a,} functional form bela actual: 
eter & COSt average 
Ry cost 


Esti- 

mated 

optimal MSE of 
replace- or 
ment 

time ~* 

1.86401 | 0.01816 | 1.17291 0.04966 
Se nee! aeee en aes 1.93160 | 0.01831 | 1.20933 0.05909 
2.525+0.06 x (i 1) 

EE ALL Cae ema es 2.01019 | 0.01897 | 1.26157 0.08345 
2.525 +0.03 x (i— 1) 

0.01918 | 1.30963 0.08491 
0.01916 | 1.41962 0.14584 
0.01997 | 1.51001 0.20026 
2.35111 | 0.02039 | 1.77728 0.37160 
ae 2.43630 | 0.02189 | 2.09759 2.63463 
1.725+0.01 x (i= 1) 

ice ie et as 2.51301 | 0.02519 | 2.60572 2.22276 
1.725+0.01 x (i= 1) 


1s | «2.58071 | 0.03489 | 3.26706 | 95.67296 
1.725+0.01 x (i - 1) 
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Figure 4. The curves of best performance {a,} sequences 











IV. ANALYSIS OF SIMULATION RESULTS 


It is clear from Tables 7, 8 and 9 that {a,} sequences with the smallest MSE when 
the shape parameter a equals 2.0, perform significantly worse when a is equal to 1.1. 
It is also clear that for a given a , the {a,} sequence performing well for small sample 
sizes will not be the best {a,} for large sample sizes, In this chapter we analyze in detail 
how the best performing {a,} sequences differ for different sample sizes, values of « for 
the underlying system life distribution, and for different costs C,, C, To make these 
detailed comparisons for i = 1, 2, ...., VY, we will examine: 


e <, the ith scheduled replacement times (censoring ages) €, averaged over 1000 re- 
petitions, 


e P(Y,<é), the proportion of unscheduled replacements from 1000 repetitions at the 
ith replacement, and 


e R,, the average actual cost per unit time for the first i replacements averaged over 
1000 repetitions. 


We begin by studying simulation results for small, moderate and large sample sizes sep- 
arately. 


A. SMALL SAMPLE SIZE N= 10 

From Table 7 on page 19, there are only three different best performing {a,} se- 
quences for V = 10. We plot the performance of all three {a,} sequences for a = 2.0, 
1,8, 1.5 and 1.3 in Figures 5, 6, 7 and 8 respectively. Figures 5(a)~8(a) are the same; 
they plot a, versus i, to give a visual representation of the three {a,} sequences. Figures 
5(b)~8(b) plot ¢, versus i, and the dotted line is the optimal age of replacement 4%. 
Figures 5(c)~8(c) plot P(X, < é,) versus i, and the dotted line is P(X, < $*), the proba- 


bility of unscheduled replacement under a policy with optimal age of replacement $*. 
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Figures 5(d)~8(d) plut 2, versus é and the dotted line is R(#*). Recall that the first 
scheduled replacement time €, = 1.0 and a, = 0.0, so that for each simulation 


, = 1.0. In Figures 5, 6, 7 and 8: 

@ the curves marked “1” 
ee 
1.400 + 0.03 x @— 1) 


e the curves marked “2” correspond to the {a,} sequence of functional form 
eae eee es 
(1.725+0.01 x (i- 1)) 


e the curves marked “3” correspond to the {a,} sequence of functional form 
et tate! eee 
1.725+0.06 x (§- 1) 


correspond to the {a,} sequence of functional form 
which performs best for a = 1.1, 1.2, ...., 1.7, 


which performs best for a = 1.8, and 


which performs best for a = 1.9, 2.0. 


1. Shape Parameter « Equals 2.0 

The average censoring time €, for all three {a,} sequences overcstimates the op- 
timal replacement time ¢* = 1.1792, for i = 2, 3, ..., 10 (Figure 5(b)). Because 
é, = 6*,+4a,, when q, is large, the scheduled replacement time for the ith replacement 
€, is also large. Thus we would expect that the {a,} sequence with larger values to yield 
larger #,, which in turn would mean that a higher proportion of replacements would be 
unscheduled. Figures 5(a), 5(b) and S(c) are consistent with this expectation. 

In Figure 5(d) R(@*) = 1.8024. This figure shows that for all three sequences 
the R, the average cost per unit time up to the ith replacement decreases with i, the 
number of replacements. This result is promising because the goal of the sequential es- 
timation procedure is to decrease costs while sampling. Even though as 1 > co, R, will 
approach the optimal replacement cost R(@*) with probability 1.0 [Ref. 8], there is no 


guarantee that R, will decrease for the first few observations. 
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Figure 5. Plot of {a,} performance when N= 10 and « = 2.0 


The ordering of the curves 1, 2 and 3 in Figure 5(d) supports the intuition that 


the {a,} sequence which yields , ‘s closest to #* will have the smallest average costs per 


unit time R. Thus, curve 3 which is closest to ¢* in Figure 5(b) is also closest to 


R(@*) in Figure 5(d). 
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2. Shape Parameter Equals 1.8 

Here ®* = 1.2032 and R(¢*) = 1.9470. Again all three {a,} sequences yield 
E> 6%, i = 2,3, n0, 10 (see Figure 6(b)). However, comparing curve 3 in Figu.2 6(b) 
with curve 3 in Figure 5(b), we see that the third {a,} sequence yields €, values closer to 
@*: for « = 1.8 than for a = 2.0. Although the third {a,} sequence clearly yields €, 
values closest to @*for « = 1.8 (curve 3, Figure 6(b)), both the second and third {a,} 
sequence have nearly identical average actual cost per unit time (curves 2 and 3, Figure 
6(d)). In fact, for « = 1.8 the second {a,} sequence performs better than the third {a,} 
sequence in that the MSE for the second {a,} sequence is smaller. 

To explain this, recall that at each replacement, we can only observe the first i 
replacements. During early stages of sampling (i.e., for small i), the information about 
the system’s life distribution will not be enough to get a very good estimates of *. 
Thus the estimator ¢* and the scheduled censoring time €, will have large variances in 
early stages of sampling. Also notice that in Figure 3 on page 15, R(«) to the left-hand 
side of the global minimum @* is steeper than on the right-hand side of *. Thus slight 
underestimates of #* can increase costs much more than corresponding overestimates 
of @%. This explains the fact that when €,>* the {a,} sequences with €, is closer to 
o*, may not do as well as {a,} sequences with slightly large é, 

We conclude that when the {a,} sequence with average censoring times 7s 
closest to @* has €,>*, and the differences the és and @* are small, then this se- 
quence might not yield the lowest average cost, and we need to check the variation in 


the ¢,’s carefully before jumping to conclusions. 
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Figure 6. Plot of {a,} performance when N= 10 and 2 = 1.8 


3. Shape Parameter a Equals 1.5 
The {a,} sequences in Figure 7 perform similarly to those in Figure 6. We see 
that the é, > o* fori = 2, 3, ...., 10 for all three sequences, and that the three curves 
of the average censoring time are even closer to #* than when « = 1,8 and « = 2.0 
(Figures 6(b) and 5(b)). The same line of reasoning used to explain Figure 6 when 


a = 1.8, can be used here to explain why the average actual costs from the three {a,} 


sequences are virtually identical (Figure 7(d)) even though the €, ’s for the three se- 
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quences are clearly separated (Figure 7(b)). In fact, the first {@,} sequence, with the 





largest €, values (curve 1, Figure 7(b)) has the smallest MSE. 
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Figure 7. Plot of {a,} performance when N= 10 and « = 1.5 


4, Shape Parameter « Equals 1.3 
Figure 8(b) shows that for all three {a,} sequences é, < o*,i = 1,2,...., 10 and 
that the first {a,} sequence (curve 1) yields &, values closest to *. Again because R(2) 
to the left of @* is steeper than R(r) to the right (Figure 3, curve 2), when several {a,} 


sequence yield €, < o* fori = 1, 2, ..., 10, then the {a,} sequence with @/s closest to 6* 
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should yicld the lowest average actual costs per unit time 2,. Thus it is not surprising 


that the first {a,} sequence has the smallest R, i = 1, 2, ..., 10 and the smallest MSE. 
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Figure 8. Plot of {a,} performance wher: N= 10 and z = 1.3 


B. MODERATE SAMPLE SIZE N=50 
We use the same approach to analyze the performance of {a,} sequence with mod- 


erate sample size as we used for small sample size. In figures 9, 10 and 11: 


e the curves marked “1” correspond to the {a,} sequence with functional form 
1 
i 
1.400 +0.03 x (f— 1) 
e the curves marked “2” correspond to the {a} sequence with functional form 
I 
(8.00 + 0.03 x (f— 1) 5 


which performs best fora = 1.1, 1.2, ...., 1.5, 


which performs best for a = 1.6, and 











e the curves marked “3” correspond to the {a,} sequence with functional form 
l 
2.225+0.06 x (i— 1) 


which performs best for a = 2.0. 
1. Shape Parameter « Equals 2.0 
Again in Figures 9(a), 9(b) and 9(c) we see the same trends that were evident for 
N = 10 when a = 2.0 (Figures 5(a), 5(b) and 5(c) on page 25). For instance, curve | 


in Figure 9(a) has the largest a, values, i = 1, 2,...., 50, in Figure 9(b) it has the largest 


Z, values and in Figure 9(c) it has the largest proportion of unscheduled replacements. 





Curve 3 has the lowest average cost and it has the smallest MSE. Thus the {a,} se- 


| 


quence of functional form ¥5354006x 0 1) 


, has the best performance for .V = 50 


and a = 2.0. 
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Figure 9. Plot of {a,} performance when N= 50 and z = 2.0 
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2. Shape Parameter a Equals 1.6 
From Figure 10(d), we sce the same phenomena that occurs for V = 10 with 
« = 1.8 (Figure 6 on page 27), ie., that although the average scheduled replacement 
times €, for curve 3 are closer to ¢* than for curve 2, curve 2 has the lower average 
cost. Again this can be explained by the fact that when é, > * the variance of the 
és has a more dramatic effect on costs when @, is close to @* than when é, is further 
away. It is interesting to note that this effect of the variation of &,’s on cost shows up 


to same extent in Figure 9 when « = 2.0. 
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Figure 10. Plot of {a,} performance when N=50 and « = 1.6 
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3. Shape Parameter « Equals 1.3 
All three {a,} sequences yield €, < b* fori = 1, 2, ., 50 (Figure 11({b)) con- 
sistent with the observations made when VY = 10 and « = 1,3 (Figure 8), the first {a,} 
sequence (curve 1 in Figure 11(b)) which yields @, values closest to @*, also has the 


smallest R,’s (curve 1 in Figure 11(d)) and has smallest MSE. 
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Figure 11. Plot of {a,} performance when N= 50 and « = 1.3 








C. LARGE SAMPLE SIZE N= 250 





Figures 12, 13 and 14 show the performance of {a,} sequences as follows: 


e the curves marked “1” correspond to the {a,} sequence of functional form 


, T.725+0.01 x(i— 1) whici: performs best for a = 1.3, 


e the curves marked “2” correspond cto the ja} sequence of functional form 
Tic.00+0.09xG= I> which performs best for « = 1.6, and 


e the curves marked “3” correspond to the {a} sequence of functional form 
3535+006x(i2 which performs best for x = 2.0. 


1. Shape Parameter « Equals 2.0 





As for the other sampfe sizes when @ = 2.0, all three {a,} sequences yield 

é > b* for i= 2, 3,...., 250. In early stages of sampling, the average costs for the second 

{a,} sequence (curve 2, Figure 12(d)) are lower than the average costs for the third {a,} 

: sequence (curve 3, Figure 12(d)). But after about 60 replacements the situation reverses 
and curve 3 ralls be’ w curve 2? in Figure 12d). From Figure 12(b), we observe that 

the average censoring times of curve 3 are large when / is small. After about 349 re- 

placements, the average censoring times of curve 3 are the closest to @*. Because the 

variance of the estimators of ¢* should be smailei, after moderate number of 1eplace- 

ments, é, of curve 3 (Figure 12(d)) decreases faster than the other two curves and fin- 


ishes lower than curve 2 (Figure 12(d)) in the end. 
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Figure 12. Plot of {a,} performance when N= 250 and « = 2.0 


2. Shape Parameter « Equals 1.6 
Although the second and third {a,} sequences are fairly different (curves 2, 3 in 
Figure 13(a)), their average costs are about the same. Before these two {a,} sequences 
intersect, the proportion of uncensored observations (i.e., unscheduled replacements) is 
higher for the second {a,} sequence than the third (curves 2 and 3, Figure 13(c)). After 


intersection the second {a,} sequence has a higher proportion of unscheduled replace- 


ment. Thus after a moderate number of replacements, it would seem that about the ‘ 


same amount of information about the system’s life distributicn and is gathered using 








either {a,} sequence, so that we could expect the average costs to be about the same. 
However, after a large number of replacements, the variance of censoring times become 


small and the {a,} sequence vielding @, ‘s the closest to the optimal replacement time 





(curve 3) has the lowest average costs. 

Although the third {a,} sequence (curve 3 in figure 12(d)) gives the lowest. Riso, 
our simulation results reveal that the variance of actual costs per unit time are larger 
than for the second {a,} sequence. Our policy is to choose the sequence with the 
smallest MSE as the best performer. Thus we conclude that the second {a,} sequence 


(curve 2), is the best performer. 
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Figure 13. Plot of {a,} performance when N= 250 and « = 1.6 














3. Shape Parameter « Equals 1.3 
Since the second and third {a,} sequences (curves 2 and 3 in Figure 14(b)) yield 
&, < * fori = 1, 2, ..., 250, we expected the sequences with €, closer to @* to have 
the lowest cost. This can be seen with curves 2 and 3, we see in Figure I4(b). The first 
sequence (Curve | in Figure 14(b)) yields &’s closest to @*. Thus it is not surprising 
to find out that the first sequence (curve 1 in Figure I[4(d)) has the lowest average cost 


and the smallest MSE. 
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Figure id. Piot of {a,} performance when N= 250 and « = 1.3 




















D. SAMPLE SIZE CHANGES 

In the last section for each fixed sample size small, moderate and large, we notice 
that as the shape parameter « decreases from 2.0 to 1.1, the best performing {,} se- 
quences tend to have larger values. This means that when the system’s life distribution 
is close to the exponential distribution, the seq. .ntial estimation procedure needs to use 
{a,} sequences with larger values to get smaller actual costs per unit time. Now we ex- 
amine how the {a,} sequences perform as sample size varies for each fix « We picked 
some of the best performing {a,} sequences from Tables 7, 8 and 9, and plotted these {a,} 
sequences in Figure 15. The short curves are the best performing {a,} sequences for 
small sample size .V = 10, the moderate curves are the best performing {a,} sequences 
for moderate sample size .V = 50 and the long curves are the best performing (a@,} se- 


quences for large sample size V = 250. 
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Figure 15, Comparison of best perforniiig (a,; when a is fixed 
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In the bottom plot of Figure 15, the {a,} sequence were the same for « = 1.4, 1.3, 
1,2 and 1.1. We already know that as the system’s life distribution approaches the ex- 
ponential distribution (i.e., as the shape parameter « decreases to 1.0) the sequential es- 
timation procedure performs best with {a,} sequences with larger values. Thus the 
condition that four different « values have the same best performing sequences, is due 
to the fact that {a,} sequences with large enough values were not tried in the simulation. 

From Figure 15, for each a, the best {a,} sequences for V = 10 are larger than the 
best sequences for V = 250. However best {a,} sequences for NV = 50 sometimes 
intersect with the best {a,} sequences for V = 250. Generally speaking, the relationship 
between three best {a,} sequences is not consistent for different « values. In early stages 
of sampling before much data been collected, the variation of censoring ages is very 
large. From Figure 3 on page 15, if the censoring time is less than the optimal re- 
placement time $*, then the censoring time will be on the steeper side of the curve. This 
means the censoring time will have a large expected average cost. Thus we would rather 
ave an {a,} sequence with large values so that the é, ‘s tend to be larger than $* to 
reduce the risk of higher costs when ¢,< 6*. 

When / becomes large, we have many observations with more information about the 
system’s life distribution which can be used to better estimate ¢*, thereby reducing the 
variance of the &,'s. Thus the {a,} sequences giving average censoring times €, closest to 
the optimal replacement time #* performed best. In general {a,} needs to decrease as 
n increases. However for distributions that are closer to exponential, @* is larger and 
the estimators 6% tend to be less than 6* , we need {a,} sequences with large values that 
decrease slowly. For large « values, 6X tends to overestimate the optimal replacement 
time, thus the {a,} sequence which drop to zero faster or have light tails, achieves lower 


cost, 
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E. COST RATIO CHANGES 
1. Optimal Replacement Time 
We recall that the expected long run costs per unit time R(s) is dependent on the 
unscheduled and scheduled maintenance cost, the system's life distribution and the 


scheduled censoring time & 


C,x Ft) + CxS 
R(t) = 1 Y) 2 x S(2) ; 
| S(x) dx 
0 
where C, and C, are unscheduled and scheduled replacement cost respectively. By ex- 
pressing R(z) as 
C ; 
wales) + S(T) 
RO=C, «( S-—— J, 


1 Stx)dx 
aR) 


we see that the optimal replacement time ¢* is a “inction of cost ratio a Thus for 
a given system life distribution, as long as the ratio is constant, then #* is the same no 
matter how large C, and C, are. Thus we cun study the effect of different costs on the 
sequential estimation procedure by changing the cost ratio. Table 10 shows optimal 


replacement times for several cost ratios. These were found using the same method used 


to construct Table 5 on page 16. 
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Table 10. COMPARISON OF OPTIMAL REPLACEMENT TIMES BY DIF- 
FERENT COST RATIO WITH E(Xj) = 2.0 


Optimal Replacement Time $* 


Shape 
Parameter « 


2.7862625 1.1846266 0.9107760 0.6948394 
2.9159756 1.2031918 0.8641901 0.7974927 


aC CE 
ee 


21.9620056 12.9815130 4.2228003 3.461828 





From Table 10 for each fixed ratio costs, the optimal replacement time $* in- 
creases as the shape parameter a decreases from 2.0.2 From Table 10 we also see that 


oan ae : ; 
for fixed « as the cost ratio —> increases, the optimal replacement time @™ decreases. 
2 


This means that when the unscheduled replacement cost is too large, we don’t want to 


risk unscheduled replacement. 


2 Note that i in Table 10, when C, = 8.0, the optimal replacement time ¢* value for a = 
1.8 is less than the ¢* value for « = 1.9. The same is true when C, = 10.0,« = 16and 1 4. 
This is caused by the specific series of Pseudo Random Numbers (PRN) used to simulate ¢*. 
When the seed for the PRN generator changes, the simulation results for these situations should 
change also. 
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2. Effect on Performance 
Let a = 1.8, the scheduled replacement cost C, = 1.0, and iet the unscheduled 
replacement cost C, have different values 2.0, 5.0, 8.0 and 10.0. To examine the per- 


formance of the {a,} sequence for different —— ratics, we simulated the sequential esti- 


ie 
mation procedure with each of the {a,} sequences from Table 6 on page 18. Again, the 


best performing {a,} sequences were chosen. They are summarized in Table 11. 


Table 11. BEST PERFORMING {a,} FOR SHAPE PARAMETER %*=1.8 






Unscheduled 
and 
Scheduled 
Replacement 
Cost 


{a,} Functional Form 

















I 
C,= 1.0 TRC re NCTE 1.725+0.01 x (i— 1) 


= 5.0 
= 1.0 eee) (i— 1) +0, 
rg: ORT ran 
C,= 10. ee re Ree ne 
C= 1.0 RICE RTGT) 


From Table 11, we observe that for C, = 8 and 10, the same {a,} sequences 









2.525+0.15 x (i— 1) 






ee 2 
2.525 +0.15x (i— 1) 


perform best for sample sizes V = 10, 50 and 250. This could be because the cost ra- 
cS 
tios Tia 8 and 10 are close, or it might be caused by the fact that we did not test 
2 
enough {,} sequences. 


4] 











We have also plotted the {a,} sequences from Table 11 in Figure 16, so that 
Wwe can make a comparison. By comparing Table 11 and Figure 16, we observe that for 
fixed sample size as the cost ratio + increases, the best performing {a,} sequence tends 
to become small. Thus the ith replacement’s scheduled replacement time will tend to 
be small and the proportion of unscheduled replacements will decrease. This limits the 


probability of an unscheduled failure. This conclusion is consistent with the conclusions 


of the previous section. 
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Figure 16. The {a,} sequence performances for « = 1.8 

















F. COMPARISON WITH FREES & RUPPERT’S PROCEDURE 

There have been two previous attempts to construct nonparametric sequential esti- 
mators of #*: Bather [Ref. 11: pp. 253-265] and Frees & Ruppert (Ref. 12: 
pp. 650-661]. Bather’s procedure was not fully sequential in that periodically, the sys- 
tem was allowed to operate until failure regardless of how long it took. Frees and 
Ruppert use stochastic approximation techniques to construct a fully sequential esti- 
mator of $* (i.e, €,< oo for i= 1, 2,....). In this section we compare the performance 
of the sequential procedure described in Chapter 3 with the perform..nce of Frees and 
Ruppert’s procedure. Because Frees and Ruppert do not estimate @* directly, we can 
not compare the procedures using mean squure error (MSE) of the estimator of ¢*. 
We can however compare MSE, the squared difference of actual costs per unit time from 
R(b*) averaged over the repetitions of each simulation, for both procedures. Because 
the ultimate goal is to reduce costs, this is, in fact, the more meaningful comparison. 

In their Monte-Carlo study, Frees and Ruppert took the unscheduled replacement 
cost C, = 5.0 and the scheduled replacement cost C, = 1.0. The system’s life distrib- 
ution is Weibull with shape parameter a = 2.2 and scale parameter 1 = 0.5. This dis- 
tribution has expected value 1.7712 and standard deviation 0.8499. The optimal 
replacement time @* = 0.99505 and optimal cost R(#*) = 1.904. [Ref. 12: p. 656] 

Frees and Ruppert update their estimator of @* after every two observations. They 
chose to look at the performance of their estimator after 10, 50 and 250 updates to re- 
flect small, moderate and large sample sizes respectively. Thus after 10, 50 and 250 up- 
dates, they used 20, 100 and 500 observations respectively. Table 12 on page 44 shows 
the performance of cost estimator from Frees and Ruppert’s sequential estimation pro- 


cedure. 
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In order to compare methods, we use the same simulation setting as Frees and 
Ruppert’s, C, = 5.0, C, = 1.0, the system’s life distribution is Weibull with shape pa- 


rameter « = 2.2 and scale parameter A = 0.5, and we select a {a,} sequence of the 
15 

(a + 50)” 
used by Frees and Ruppert. The performance of our cost estimator is shown in the 


Table 13. 


functional form . This {a@,} sequence roughly corresponds to the sequence 


Table 12. THE PERFORMANCE OF COST ESTIMATOR FROM FREES AND 
RUPPERT’S PROCEDURE 
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Table 13. THE PERFORMANCE OF COST ESTIMATOR FROM OUR SE- 
QUENTIAL ESTIMATION PROCEDURE 
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Recalling that Frees and Ruppert’s need two observations for each update of their 


estimator of @ and we only need one observation to update our estimator of *. The 


results of our simulation after 20, 100 and 500 observations are comparable to Frees and 


Ruppert’s simulation with 10, 50 and 250 updates of $* respectively. Comparing 


Table 13 to Table 12 with the same number of updates, the variance of our simulated 


average costs are slightly higher than Frees and Ruppert’s simulation results. Elowever: 


Our average costs are much lower than Frees and Ruppert’s average costs. 


For small sample size (20 observations), for this {a,} sequence , the MSE for both 
procedures are about the same. For large sample size the MSE of our procedure 
is smaller than that of Frees and Ruppert, even using half the number of observa- 
tions. 


This choice of {a,} gives the best performance for Frees and Ruppert’s procedure 
and was chosen knowing that the underlying distribution is Weibull (a= 2.2, 
A=0.5). It is probable that a different choice of {a,} would yield better results for 
our sequential procedure. 


The intermediate estimators of Frees aiid Ruppert’s sequential estimation proce- 
dure are complicated to calculate. 


Thus we may conclude that our sequential estimation procedure for age replacement 
policies is competitive to Frees and Ruppert’s procedure, and is in many ways better. 

















V. GRAPHICAL DETERMINATION OF THE SCHEDULED 
REPLACEMENT TIME 


A graphical method can be used to estimate the optimal age replacement policy 
based on observing a sample of iid. system lifetimes (Ref. 3: pp. 113-115]. This 
method obtains the estimate of 6* froma scaled total time on test plot. The advantage 
of this method is that you can see immediately how sensitive the estimate of 6* is to 
change in the cost ratio. This method can be extended to estimate @* based on data 
from our sequential estimation procedure. In this chapter, we describe the graphic de- 
termination method for iid. data and apply it to the 22 D9G-G66A Caterpillar Tractor 
Engines (this data was used as an example in Chapter 2). We also indicate how this 


method can be extended to our sequential estimation procedure. 


A. SCALED TOTAL TIME ON TEST PLOT 

Plotting the data is often the first step to unlocking information contained in data 
about the underlying model. For example, a plot of the empirical distribution function 
contains information about the probability density. If the empirical distribution looks 
concave, then the density may be decreasing. However, this plot is scale dependent so 
that the perceived shape of the plot depends on the choice of plotting scale. 

A total time on test plot provides information about the failure rate. Analyzing 
failure data, it is often the failure rate which is of chief interest. For example, if the 
failure rate is constant or decreasing, we know that we should not adopt a scheduled 
replacement policy, since an old system in this case is actually “better” than a new sys- 
tem. Consider an experiment where a systems system lifetimes Xy Xqy vey X_ are OD- 
served, Let n(u) be the number of systems that survive to age u, then the total tirae on 


test to age x is defined to be 
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rx 
T(x) = | n(u) du. 
0 
Let Xy < Xy Su. S Xi, be the order statistics of X,, X,, ...., X,, then 
Xo 
T(Xp) = n(u) du 
0 
= nX v1) + (n = DC ‘Q) = Xi) tou. oF (nt —i+ NX nest X-1)) 


is the total time on test to age \,. We call 


_ T(x) 
TX) 
' T(X,) oe . , 
the scaled total time on test at age x. A plot of T,) versus the empirical distrib- 
() 


ution F(X) = 7 fori = 1, 2, ...., is called the scaled total time on test plot. This plot 
provides information about the failure rate. If the total time on test plot is strongly 
concave, then there is evidence that the underlying distribution is IFR and a scheduled 


age replacement policy makes sense. 


B. OBTAINING THE AGE REPLACEMENT TIME FROM THE PLOT 

Let C, be the cost (in dollars) of an in-service tractor engine failure replacement and 
C, be the cost (in dollars) of a scheduled tractor engine replacement, where C, > C,. 
After 22 observations, ¢* can be estimated from the tractor engine failure ages using the 
procedure described in Chapter 2 with a, = €, = oo, i = I, 2,...., 22. It can also be 
estimated using the scaled total time on te . plot. Table 14 on page 48 shows the or- 
dered actual failure age of 22 D9G-66A Caterpillar Tractor Engines, and the empirical 


distribution of tractor engine failures is plotted in Figure 17 on page 49. 


(C, - G) 
the scaled total time on test plot through this point. Then drop down to the x-axis from 


To estimate $* plot the point on the x-axis and draw a tangent line to 
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Table 14. THE ORDERED AGE OF FAILURE TRACTOR ENGINES 
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: i f : 
the point of tangency closest to the value - when j, =1, 2,....,. Then X(,, will be the 
estimated optimal replacement age. A straight forward geometric argument shows that 
estimating * from the scaled total time on test plot is identical to the estimate obtained 


by minimizing R(T) [Ref. 3: p. 114]. 
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Figure 17, The empirical distribution and the total time on test plot for Caterpillar 


Tractor Engine data 


C. THE TOTAL TIME ON TEST PLOT OF TRACTOR FAILURE ENGINES 
The bottom plot in Figure 17 is the scaled total tir‘e .n test plot when the un- 
scheduled replacement cost C, = 200 and the scheduled replacement cost C, = 100. 


On the horizontal axis we can find the point = -1.0. We draw a tangent line 


we 
(C, - G) 
to the scaled total time on test plot through this point anc drop down to the x-axis from 


the point of tangency, we find the value = = 0.273 = s. Then the order statistic X% 


_(i.e., 5085 hours), is the estimated optimal replacement time. Changing the unscheduled 


replacement cost C, to 400 and 700, and holding the scheduled replacement cost C, at 
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100, from the scaled total time on test plots in Figure 18 on page 50, we see that the 
estimated optimal replacement age is 4394 hours and 2690 hours respectively. The re- 
sults obtained from this figure are consistent with the optimal replacement ages esti- 
mated using the sequential estimation procedures in Table 4 on page 11. The value of 
the estimator using the scaled total time on test plot depends upon how closely the es- 


timated distribution from the sample approximates the true underlsing life distribution. 
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Figure 18. The total time on test plot for Caterpillar Tractor Engine data 
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We note that this procedure can be adapted to the sequential estimation procedure 
(Ref. 8} by replacing the empirical distribution as an estimator of F, by the product-limit 


estimator of F described in Chapter 2 and defining 


| § (u) du 
T(x) 0 


TZ) i ‘ 








S(u) du 
0 

where 

¢ _ Pe (Sivas) Gee a) 

Sy) = (251) 

(t:Zy <u) 
. TEES: este hope ee T(Zy) " A 
is the product-limit estimator of S=1—F. Then plot T(Z,.) versus F(Z) = 1 — S(Z) 
(n) 


to get the scaled “total time on test” plot. Following the same steps outlined for the 


2 


iid, case, plot Gao E the x-axis and draw the tangent line to the scaled total time 
Pp? 








C, Fs ‘ . A . 
on test plot through the point G C . The estimator of * is F-'(x,) where x, is the 
fo 2 
value on the x-axis, below the point where the tangent line intersects the scaled total 


time on test plot. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


Under an age replacement policy, a system is replaced at failure or after being in 
service for ¢ units of time, whichever comes first. The time r is called the scheduled 
censoring time or the scheduled replacement time. An is optimal replacemeni time @* 
achieves the smallest long run expected cost. An important problem is the estimation 
of @* when the form of the system's underlying life distribution is unknown. We show 
by example and through simulation that substantial cost savings can be effected using 
the sequential estimation procedure for @* described in Chapter 2. 

An important part of our analysis is to provide practical guidelines for choosing 
é,= $*, +4, the schedv'ed replacement time to be used for the ith system, where $%, is 
the estimator of 6* based on (i— 1) replacements, and {a,} is a fixed sequence of con- 
stants. From our analysis, in the early stages of sampling we would rather choose large 
values for a, so that &, tends to be greater than @*. If the scheduled replacement times 
€, are too small in the early stages, then most of the system life times are censored. Thus, 
very little information is collected about the underlying distribution. This leads to poor 
estimates for @*. After a large number of replacements, the estimates of the underlying 
distribution and therefore * are relatively good. A good age replacement policy will 
keep the scheduled censoring time closely approximate the optimal replacement time 
p*. Therefore the best performing {a,} sequences have smaller a, for large i. 

The proper choice of {a,} sequence is a very difficult. The best choice of an {a,} 
sequence for sequential estimation of age replacement policy depends on the system’s 
underlying life distribution which is generally unknown, the costs C,, C,, and the sample 
size .V at which sampling stops. From our simulation, we see that when the underlying 


distribution is close to exponential the {a,} sequences should have larger values and 
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should not decrease as fast as those whose failure rate increases more rapidly. For large 
sample sizes {a,} sequences should decrease fuster in early stages than small sample sizes. 
In addition as the cost ratio 4 increases, the best performing {a,} sequence have 
smaller values. 


Future research might be concerned with the following subjects: 


¢ Sequential estimation \ hen the underlying life distrivution F comes from a 
parametric family of distributions. 


¢ How to use the data to estimate the {a,} sequence. 


e If an imperfect repair model is permitted, then the next failure depends on system 
history, failure distributions are not independently and identically distributed 
(i.i.d.), how to change our Sequential Estimation Procedure to fit new situation. 


¢ Minimizing long run expected costs per unit time gives the same optimal ¢*for 
C, = 2, C, = | and costs for C, = 2,000,000, C, = 1,000,000. Clearly minimizing 
long run expected costs is not appropriate under all circumstances. Other cost 
functions need to be considered. 
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